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Abstract 

The evolution of scalar, electromagnetic and gravitational fields around spherically symmetric 
black hole surrounded by quintessence are studied with special interest on the late-time behavior. 
In the ring down stage of evolution, we find in the evolution picture that the fields decay more 
slowly due to the presence of quintessence. As the quintessence parameter e decreases, the decay of 
Q ■ £ = mode of scalar field gives up the power-law form of decay and relaxes to a constant residual 

^ • field at asymptotically late times. The £ > modes of scalar, electromagnetic and gravitational 

fields show a power-law decay for large values of e, but for smaller values of e they give way to an 
exponential decay. 
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I. INTRODUCTION 



Black holes are among the simplest objects in the Universe as they can be fully described 
by merely three quantities viz., mass, charge and angular momentum. During the process 
of gravitational collapse of a charged rotating star leading to the formation of a black hole, 
the geometry outside the star relaxes to the stationary Kerr-Newman spacetime. This 
result, what is now known as ^^no-hair theorem!^ means that any "hair" other than its mass, 
charge and angular momentum, will disappear after the collapsing body settles down to its 
stationary state. The mechanism responsible for the vanishing of these quantities remained 
unknown till Richard PriceP] who, making a perturbative analysis of the collapse of a 
nearly spherical star, showed that for a field with spin, s, any radiative multipole {I > s) 
gets radiated away completely, in the late stage of collapse. Further, he showed that at late 
times the field dies out with a power-law tail where p = 1 if the multipoles were 

initially static and p = 2 otherwise. Later, in[2| a detailed analytical analysis supported 
by numerical results of massless fields around Schwarzschild and Reissner-Nordstrom(RN) 
black hole showed that power-law tails are also present on null infinity where they decay as 
n-i^+p) g^Yid on future event horizon as 

The studies of the evolution of charged scalar field around RN spacetime reveal that a 
charged hair sheds slower than a neutral hair. In contrast with the massless fields, massive 
fields with mass, m, have an oscillatory inverse power-law behavior t~^^^^^'^^ sin(mt), at 
intermediate late times |5[. But it was shown that in the asyniptotic late times another 
pattern of oscillatory tail of the form, t^''^/^^ sin{mt) dominatesjo, 7|. The late-time tails 
appearing in the propagation of fields in the higher dimensional black hole were studied in 



a- 

The evolution of fields propagating on spacetime with a non vanishing cosmological con- 



stant were addressed in 
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3l-|l4j and they have demonstrated the existence of exponentially 
decaying tails at late times contrasting the power-law tails in asymptotically flat situation. 

Earlier studies provide us the following picture of evolution of perturbations in the black 
hole spacetimes. It involves three stages 15|], the first one is an initial response, determined 
by the particular form of the original wave field followed by a region dominated by damped 
oscillation of the field called quasinormal modes(QNMs), which depend entirely on the back- 
ground black hole spacetimes. The QNMs play an important role in almost all astrophysical 
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processes involving black holes and the spectra of QNMs are expected to be detected by the 
future gravitational wave detectors. At late-times the QNMs are suppressed by the so called 
tail form of decay. This comprises the last stage of evolution and is an interesting topic of 
study, since it reveals the actual physical mechanism by which a perturbed black hole sheds 
its hairs. 



In the past two decades there have been growing observational evidences [ly, [iTj which 
indicate clearly that our universe is expanding in an accelerated pace, indicating the presence 
of some mysterious form of repulsive energy called dark energy. In order to explain the nature 
of dark energy, several models have been proposed. The simplest option being Einstein's 
cosmological constant. A, which has a constant equation of state with state parameter, 
e = — 1, but it needs extreme fine tuning to account for the observations [3]. Models were 
proposed, replacing A with a dynamical, time-dependent and spatially inhomogeneous scalar 
field now called quintessence, which can have an equation of state, — 1 < e < — l/sjiol, Q. 

Spherically symmetric black hole solutions 2l|] were obtained in quintessence model of 
dark energy and many authors 22|-[30|] considered the evolution of various fields around black 
holes surrounded by quintessence. All these studies are based on the calculation of the QNMs 
using the third order WKB approximation method while the late-time behavior remained 
unexplored till now. In fact the late-time decay is determined by the asymptotic curvature 
of the spactime, so it is interesting to see how the fields evolve in a spacetime in which the 
asymptotic structure is determined by the quintessence field and our study fills this gap. 

The rest of the paper is organized as follows. In SectHIlwe introduce the master wave 
equation for massless scalar, electromagnetic and gravitational field perturbations around 
black hole surrounded by quintessence. The numerical method used to study the time 
evolution is explained in Sect JIIII and the results are presented. The conclusions are given 
in SectHVl 



II. FIELDS AROUND BLACK HOLE SURROUNDED BY QUINTESSENCE 

The exact solution of Einstein's equation for a static spherically symmetric black hole 
surrounded by the quintessential matter, under the condition of additivity and linearity in 
energy momentum tensor, was found in|2l|. 
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^ -f(r)df + f{r)-V + r^de'' + sin^ ^^0"), (1) 

where /(r) = (^1 — ^ — ^3^+1 ^, M is the mass of the black hole, e is the quintessential state 
parameter and c is the normalization factor which depends on the density of quintessence as, 
Pg — ^ ^3(i-e) ■ It is difficult to analyze the above metric for arbitrary values of the parameter 
e. For our study we take three special cases of the quintessence parameter, e = —1/3, —2/3 
and —1 along with the Schwarzschild(Schw) case(c = 0), so that the calculations become 
viable. 

Casel e = —1/3 : the black hole event horizon is located at r^i = 2M/(1 — c), corre- 
sponding to zero of the function, /(r) and the tortoise coordinate can be defined as. 



r* = r + reiln{r - rgi). 



(2) 



Case2 e = —2/3 : here, in addition to the black hole event horizon at r = re2 the 

spacetime possesses a cosmological horizon at r = rc2, with < rc2- In terms of these 
horizons we can write /(r) as. 



/(r) = -(r - re2){rc2 - r). 



(3) 



The surface gravity associated with the horizons at r = rj, is defined by 

1, 



i^i = 2\^f/dr\ 



and we get. 



«e2 



c{rc2 - re2) 



2r, 



e2 



These quantities allow us to write. 



«c2 



c(re2 - rc2) 



2r, 



(4) 



(5) 



c2 



+ 



/ 2Ke2{r - re2) '2Kc2{r - rc2) ' 

Now we can get an expression for the tortoise coordinate, = J f'^dr as. 
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CaseS e — —l: this is the extreme case of quintessence, the Schwarzschild-de Sitter (SdS) 
spacetime. The roots of the polynomial equation, /(r) = corresponding to the event 
horizon at r = Tgs, the cosmological horizon at r = and a negative root at r = tq = 
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~{i^e3 + fcs), with ro < < Tcs- We can define the tortoise relation in terms of these 
horizons and the corresponding surface gravity as|13|, 
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(8) 



The evolution of massless scalar field $, electromagnetic(EM) field and gravitational 
perturbations(GR) h,,,,^ in the spacetime qy„y, specified by Eq.([T]) are governed by the Klein- 



Gordon 



22| . Maxwell [311 and Regge- Wheeler |24j. |32| equations, respectively, 



F.r = 0, 



with F^jj = 



A —A 



(9) 
(10) 

(11) 



where R^uid + h) is the Ricci tensor calculated by considering small perturbations, h^u, in 
the background spacetime g^^. For gravitational field, there are two modes of perturbations, 
axial and polar. Here we are considering only the axial perturbation for our study. The two 



were shown to have a similar QNM spectra and late time behavior in de Sitter spacetime 
The radial part of all the above perturbation equations can be decoupled from their angular 
parts and reduced to the form. 



-Q^+g^)Mt,r) = -Veir)Mt,r) 



(12) 



where the tortoise coordinate is defined by, dr^ 
given by. 



^dr and the effective potentials are 



/ 



Vsc 
Vgr 



fir) 



+ _ 6M _ c(3e + 3)~ 

3e+3 



(13) 
(14) 
(15) 



III. NUMERICAL INTEGRATION AND RESULTS 



The complex nature of the potentials makes it difficult to obtain the exact solutions of 
Eq. (fT2l) and we have to tackle the problem by numerical methods. A simple and efficient 



method to study the evolution of field were developed in[2], after recasting the wave equation, 
Eq.f lT2|) . in the null coordinates, u = t — r^. and v = t + as. 



dudv 

Using the following discretization. 



Tp{u,v) = V{u,v)ijj{u,v). 



(16) 



(17) 



we perform the numerical integration on an uniformly spaced grid with points, N{u + 
/i, v+h), W{u+h, v), E{u, v+h) and S{u, v) forming a null rectangle with an overall grid scale 
factor of h. Since the late-time behavior of the wave function is found to be insensitive to the 
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initial data, we set ipi.'^^''^ = 0) = and a Gaussian profile, ipi."^ = 0, t;) = Aexp 
In all our calculations we set the initial Gaussian with width, a = 3 centered at vq = 10. 
Due to the linearity of equation one has the freedom to choose the amplitude of the initial 
wave A = 1. To proceed the integration on the above numerical scheme one has to find the 
value of the potential at r(r,,) = r{{v — u)/2) at each step. Employing the Newton-Raphson 
method one can invert Eqs.([2]), ([7]) and (|8]) and get r(r^,) for the characteristic integration. 
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FIG. 1: Log-log graph of the evolution of scalar field in a quintessence filled black hole spacetime 
with quintessence parameters e = —1/3, —2/3, —1 and c = 10^^/2, in comparison with that in the 
pure Schwarzschild spacetime. The case of £ = mode is shown in the left plot and £ = 1 on the 
right. 



Figure [T] shows typical evolution profile of scalar field in a quintessence surrounded black 
hole, in comparison with that in the pure Schwarzschild spacetime. We observe that the 
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evolution shows deviations from the Schwarzschild case, after initial transient phase. The 
damped single frequency oscillation(QNM) phase and the late-time tail of decay in the final 
phase shows the characteristics of the quintessence. We analyze each phase in detail. 



A. Quasinormal modes 



QNMs of various field perturbations around black hole surrounded by quintessence were 



evaluated in 



22 



30| using third order WKB method. We are not going for a quantitative 



study of QNMs but the typical nature of QNM phase can be seen in Figure [21 where we 
have plotted the evolution profiles of the different fields in a quintessence surrounded black 
hole for £ = 2, in comparison with the corresponding mode in the pure Schwarzschild black 
hole. On the left pannel the decay of scalar field is shown for e = —2/3 and —1 along with 
the Schwarzschild case. 




50 100 150 200 250 300 



FIG. 2: Semi-log graph of the QNM phase of evolution for the i = 2 modes of various fields. The 
plot on the left is for scalar field where the curves from bottom to top is for the Schwarzschild, 
e = —2/3 and -1 cases, respectively. On the right EM and GR fields are shown where bottom curve 
in each sub plot represents the Schwarzschild case and the other ones are for e = —2/3. We take 
c = 10^/2 for these plots. All the three fields damp more slowly in the presence of quintessence. 



It can be seen that the decay of the fields in the QNM phase is slower in the pressence 
of quintessence than in the pure Schwarzschild cases. Also the damping is lower for smaller 
values of e. This observation is in contrast with the results presented in 22|], where thay 
obtained that a massless scalar field damps more rapidly due to quintessence. Electromag- 
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netic and gravitational field perturbations also show similar behavior as that of scalar field, 
in the QNM phase. So we have shown the e = —2/3 case only in the plots. These results 
are in agrement with the observations made in 23|-|25[|. using the WKB method. 



B. Late-time tails 

The QNM phase is followed by the regime of late-time tails of field decay. We find that 
the nature of late-time tails are sensitive to the parameters of quintessence, the angular 
momentum index and also the type of field perturbations under consideration. Below we 
discuss the different cases in detail. 

1. ^ = mode 

Scalar field perturbations can have the i = mode, while EM and GR fields can only 
have modes with i > 1 and i >2 respectively. From Figured! which is a log-log plot, it can 
be seen that for the i = mode, scalar field with e = — 1/3 case of quintessence has the form 
of a power-law tail, with slightly slower decay rate than the corresponding Schwarzschild 
tail. 

For e = —1/3, the effective potential of scalar field is positive definite for G [— oo, +oo] 
and have a potential barrier near the event horizon but vanish asymptotically as — )■ ±oo. 
But for smaller values of the quintessence parameter, e = —2/3 and -1, after a barrier nature 
near the event horizon, the potentials with i = mode of scalar field, vanish at some r* = r° 
and there after form a negative well in the range r° < r,, < +oo. Figure [3] shows potentials 
experienced by the scalar, EM and GR fields. To clearly see the behavior of the potential 
between the horizons, we choose the value of the parameter c = 10~^ for e = — 1 case and 
a larger value of c = 10^^/2 for e = —2/3 case, so that we can bring down the seperation 
between the horizons. 

This peculiar shape of potentials for e = —2/3 and -1 reflects in the late-time tails, 
that the fields settle to a residual constant value, rather than decaying to zero. We have 
evaluated the field, = ip/r on the surface of constant r(at r^, = 10), the cosmologi- 
cal horizon(approximated by the null surface v = Vmax) and the black hole event hori- 
zon(approximated by the null surface u = Umax)- In the asymptotic late times field settles 
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FIG. 3: Plot of effective potentials 
for e = -2/3 with c = 10-2/2 and 
e = — 1 with c = IQ-^. Curves from 
bottom to top is for £ = 0, 1, 2 and 3 
modes. The scalar (solid curves), EM 
(dashed) and GR(dotted) fields are 
drawn. At large distances, GR po- 
tential is merged with the scalar case 
for e = —2/3 and with EM case for 
e = —1. The potential is scaled as 
V* = V{re-rf{t + l/2), to enhance 
the nature at large r. 
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FIG. 4: The decay of £ = 
mode of scalar field with 
quintessence parameters e = 
-2/3 with c = 10-V2, and 
e = -1 with 10"^ The fields 
are evaluated on the constant 
surface = 10(dotted curves), 
black hole event horizon(dashed 
curves) and cosmological hori- 
zon (solid curves). 



to the same constant value, (^o on all these surfaces, as shown in Figure HI It was first 
observed in[l^, for scalar field in the SdS spacetime and they have shown that at late times, 



0|,=o^0o + 0i(r)e-2'^^*. (18) 

Figure El demonstrates the £ = mode of scalar field for different values of c. The smaller 
the value of c, the later the field descends to the constant value, 0o- The dependence 0o, on 
the parameter c is shown in Figure |6l We observe that the asymptotic value of field scales 
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as ipo ~ c^-^'^^, when e = —2/3 and ipo ~ c"'^^^, when e = — 1 as found in 12|. The possible 
numerical errors in these results can be visualized from the right plot in Figure |6l The error 
increases quickly as the grid scale, h increases. 
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FIG. 5: The decay of £ = mode of scalar field. On the left pannel, the e = —2/3 case is shown 
where the curves from top to bottom represents the plots for c = 10^^/2, 10^^, 10~^/2 and 10^'^, 
respectively. The e = — 1 case is shown on the right, where the curves from top to bottom represents 
the plots for c = 10-^/2, IQ-^, 10"'*, 10"^ 10"^ and 10"'^, respectively 
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FIG. 6: On the left pannel, the asymptotic value of scalar field, <j)o versus c is plotted in logarithmic 
scale. A linear fit(dotted line) gives the slopes 1.873, for e = —2/3 and 0.995, for e = —1. On the 
right ^n|(/>o| for different numerical grid resolution is plotted for the e = —1 case. 
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2. £ > modes 

Now we report the results obtained for the decay of ^ > modes of field perturba- 
tions. Figure [7] demonstrates the evolution of scalar and electromagnetic and gravitational 
fields around a black hole surrounded by quintessence with e = —1/3, along with pure 
Schwarzschild case. The fields, are evaluated on the surface r^, = 10. A straight line 
in such a plot shows a power-law tail. The late-time tails follow the power-law decay for 
e = —1/3, but with a smaller decay rate, than the Schwarzschild case of, ■?/' ~ t"*^^^"'"'^^ For 
the quintessence case with c = 10~^/2, we get i\) ~ ^-(2^+2-7) _ 

FIG. 7: Log-log plot of the evolu- 
tion of fields for e = —1/3, with c = 
10^^/2. On the left decay of scalar 
field is shown in comparison with the 
Schwarzschild case. Curves from top 
to bottom is for i = 0,1,2,3 and 
4. The decay has a lower power-law 
constant, t~(2^+2.7) £qj. quintessence 
case. The i = 3 mode of different 
fields is shown on the right. All the 
field has an identical tail of decay for 
e = -1/3. 

As the value of the quintessence parameter decreases to e = —2/3 and -1, though in 
the intermediate time the fields follow the power-law decay, in the asymptotic late-time 
it deviates from the power-law form. Figure [S] is a semi-log plot of the decay of scalar, 
electromagnetic and gravitational perturbations for the e = —2/3 case with c = 10"^/2. A 
straight line in such a plot corresponds to an exponential decay. We can see from the plot 
that for i > 0, the intermediate time power-law form of decay is replaced by an exponential 
decay of field, in the asymptotic late-time. In particular, in the asymptotic late times, the 
decay can be fitted in the forms. 
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■0^6 P^''^^\ for Scalar and GR, 

where p is some constant, for the plots shown in Figure [H] with c = 10~^/2 we get 
p = 1.089. The evolution of field for the e = — 1 case is plotted in Figure [9l in semi-log scale, 
where we have chosen c = 10~^, in order to see the decay for large i. All the £ > modes 
at late-time relax as an exponential tail, according to the form, 

iP ~ e"^"^^*, for Scalar, 

^ ^ g-fcc3(^+l)t^ f^^ EM (20) 
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FIG. 8: Semi-log plot of the decay of 
scalar and electromagnetic and grav- 
itational perturbations for e = —2/3 
with c = 10^^/2. GR has an identical 
decay pattern as that of scalar field. 



The late-time tails in black hole spacetimes are understood as a result of particular form 
the effective potential far away from the black hole. The tails are generated by the scattering 
of the waves off by the effective potential at great distances. Performing a systematic study 
of various potentials Ching et al. 33| and later by Hod 3J] had given the following heuristic 
picture. The late-time tails reported by an observer at r^, are caused by the waves that are 
scattered by the potential, V{r'J at the point r'^ ^ r* and at late times, ip V{r'J ~ V{t/2). 
Their studies show that the inverse power-law decay is a property of asymptotically flat 
spacetimes, and that different behaviors should be expected in black-hole spacetimes that 
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/ merged with the EM profiles 

have different asymptotic properties. The results that we obtained can be explained in this 
view. For e = —1/3, all the fields have inverse power nature of the effective potentials as 
that in asymptotically flat spacetime and lead to inverse power-law decay of fields at late 
times. When e = —2/3 and -1 the potentials exponentially drop off in the asymptotic region, 
V{r^) ~ Q-^c^* and one can expect an exponential tail of the forms Eg. (1191) and (!20|) at late 
times. 



IV. CONCLUDING REMARKS 

Assuming a quintessence model of dark energy as the cause for the observed accelerated 
expansion of the universe, the paper studies the evolution of various field perturbations, 
especially the late-time decay of wave fields around a black hole spacetime. We find in the 
evolution picture that in the QNM stage, all the fields decay more slowly due to the presence 
of quintessence. At late times, QNMs are suppressed by the tail form of field decay. As 
the value of the quintessential parameter e, decreases, the late-time decay of £ = mode of 
scalar field gives up the power-law form of decay, relaxing to a constant residual field. This 
asymptotic value of the field is determined by the value of the parameters, e and c. 

For the behavior of £ = mode when e = —2/3 and -1, there is no analogue case in 
black holes with asymptotically fiat spacetime. This can be understood as a consequence 



of cosmological no-hair theorem ll[|. Comparing the situation in SdS spacetime with the 



scattering inside the charged black hole it was argued that a constant mode can be trans- 
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mitted to both black hole and cosmological horizons. However the constant solution carries 



12| 



no stress-energy. A theoretical justification for the behavior of £ = mode is given in 
in SdS spacetime. 

For large values of e, the C. > s modes of scalar, electromagnetic and gravitational per- 
turbations, still show a power-law decay, having a slower decay rate than the corresponding 
Schwarzschild case. As the value of e decreases, the power-law decay gives way to an expo- 



nential decay. The results for the SdS black ho^ 



are consistent with the previous studies |12|-|l4|. 



e spacetime, the extreme case of quintessence 



Our results may give insights to the issues related to the instability of Cauchy horizon of 
the black hole. For instance, it was demonstrated in [35] that for all values of the physical 
parameters, the Cauchy horizon inside charged black holes embedded in de Sitter spacetime 
is unstable to linear perturbations. In view of the above results, it deserves a further 
detailed study on the radiative tails of perturbations around the charged black hole in 
the quintessence model which is under investigation. 
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